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II 

)Xhi8 study relates an attempt to match the HCMM- & U2HCMR>derived 
temperature data over two testsltes ( of very local size) to similar data 
collected in the field at nearly the same times. Considerable logistical 
problems were encountered. The results indicate that HCMM investigations 
(using resolution cells of bOOm ,or so) are best conducted with areally^ 
extensive sitae, r/ither than point observations. The DAY-VIS imagery is of 
of e: cellent quality, ai.d has considerable usefullness for GEOLOGY , especially 
for structural (lineament) studies. For thes® purposes one does not need 
the Day-Night registered Imagery, except that as a single product, not to be 
used for further calculations, the iCELTA-T Imagery is most useful, again 
for structural geology. Our attempts to register the ground-observed temps. 

I (even for 0.5 sq.mile targets) were unsuccessful, due to the excessive 
pixel-to-pixe.'l noise on the HCMM data. 

(2) Several covipurer models were explored, and related to changing of^ 
the values of-t>te?sial parameters, with observed data. Unless quite complex 
models, with many parameters which can only be observed (werhaps not even 
measured!) under remote sensing conditions (e.g. roughness , wind shear etc.) 


be most 

II. Kar mm nmmmt mr mmrnUH 

THERMAL INERTIA, TEMPERATURE, 

DIURNAL TEMPERATURE, HCMM. 

14 OMkuMK immmm 

11. Umttf 

llnelatilfled 

» OmmtmOmt.i/mommoh 

IfcteURtlfiad, — 

91. Na.al Sipa 

144 

n ^ 



2 


TABLE OF CONTENTS 

P£ae 

ABSTRACT 1 

I. INTRODUCTION 5 

II. INITIAL FIELD STUDY - YERINCTON NEVADA (AUGUST 8-9, 1977). . 6 

A. Field Test Sites 7 

1. An&conda Dump Site Station 7 

2. HacArthur Station 8 

6. Results 6 

1. Soil Temperature Measurements 8 

2. Sol-a-Meter Readings 22 

3. Net Radiometer 

k. Exotech Landsat Band Radiometers 26 

5. Recording Thermometers 27 

6. Soil Samples 27 

III. THERMAL PARAMETER DETERMINATIONS AND MODELLING (AUGUST 

MISSION) 32 

A. Thermal Model Parameters 33 

B. Test Site Model Results 35 

C. Approximation to a Radiant Heating Method for 

Determining Thermal Inert’.? 43 

IV. SECONDARY FIELD STUDY - YERINGTON NEVADA (DECEMBER 8-9, 

1978 44 

A. Description of Test Siic; 44 

B. Thermal Data 45 

C. Model Results. . 53 


3 


TABLE OF CONTENTS (Continued) 


Pioe 


V. P-3 (MMS), U-2, AND HCMM DATA INTERPRETATIONS 57 

VI. CONCLUSIONS 74 

VII. REFERENCES 75 

VIII. APPENDICES (See Part 2 of 2) 

1. Data Sheets 76 

2. Calibration of the Sol-a-Meter 92 

3. Recording Thermometers 96 

h. Soli Samples 109 

5. Calibration of the PRT-A and PRT-5 units 113 

6. Photographs of the Field W&rk 126 


7. Report on Thermal-Inertia Related Computer Programs 

FS6 & FSFL2 (CSIROTEMP) .WATSON and SURTEMP 131 





5 


I. INTRODUCTION 

The understanding and interpretation of thermal infrared data has 
continued to be a significant goal of several research organizations 
throughout the worlo. The Stanford Remote Sensing Laboratory, one of 
these groups, is keenly aware of the difficulties that still exist in the 
determination of thermal parameters at the ground>air interface. The re* 
lationship of these parameters to thermal mapping employing techniques 
was a major aspect of this study. 

The work performed under this contract attempts to assess the type and 
complexity of the thermal models required for rock and soil parameter dis* 
crimination employing aircraft and satellite (HCMM) thermal infrared data. 

A field me(asurement program was performed to compare the modelling results 
and assess their accuracy. This field data was then used to establish 
local calibration sites to which tiie aircraft and satellite data could be 
related. 

The ultimate objective of the study was to establish the feasibility 
and value of thermal infrared data in the delineation of hydrothermal I y 
altered areas in the Yerington Nevada test site. The spectral filtering 
of the current Landsat system does not permit the separation of ferric 
Iron hydrothermal alteration from ferric iron coatings on unaltered vol- 
canic rocks. However, the marked density differences between the hydro- 
thermal I y altered and unaltered rocks should produce, after the proper 
modelling of surface me. logical variables, a resolvable difference 
in the thermal property of these units. The study, therefore investi- 
gated the combined value of Landsat and HCMM data in the discrimination 
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of hydrothermo) alteration zonef at Yerington Nevada. 


II. INITIAL FIELD STUDY - YERINGTON NEVADA 
(AUGUST 8-9, 1977) 

2 

Field measurements coincident with U-2 and M S data acquisition were 
made August 8-9t 1977 at two sites near Yerington, Nevada. The mission was 
designed to Investigate the relationship of thermal parameters to rock 
density for mineral exploration, and to further test the validity of ther- 
mal n»delling by relating overflight data over Yerington to the two local 
cal ibration sites. 

In order to meet these objectives, many types of data were collected 
every 38 minutes over a 2k hour period. Temperature of the soil surface 
was recorded using thermocouples, and PRT-A and PRT-5 radiation thermome- 
ters. The field team also recorded soil temperatures at various depths, 
net radi.>tion, and short -wave length radiant flux incident upon the ground 
surface. Exotech, (Landsat band), radiometers were employed to measure 
surface albedo. Because the thermal properties. Inertia (K^)^ and dlf- 
fusivlty (K/|^c), of the surface materials are very dependent on the local 
meteorology, several additional measurements were made. Data recorded to 
describe the local meteorological conditions Included the air temperature 
just above the soil-air interface and approximately one meier above the 
soil, the near»surface humidity, wind velocity, and the percent cioud 
cover. Soii moisture samples were taken at depths corresponding to the 
soil temperature probes. 

Additional calibration measurements were made at a leach pond near 
the Anaconda Dump Station. Two recording thermometers were used to record 
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the subsurface water temperature continuously during the 2k hour data col* 
lection period. 

A. FIELD TEST SITES 

1 . Anaconda Dump S I te 

The main data recording station was located on the Anaconda Company's 
waste rock dump north of the Yerington open^plt mine near a level access 
road (Plate 1). The site is bordered by waste rock mounds about 8 feet 
high, and on the south side of the site there Is a cliff facing the leach 
ponds to the north (plates 2 and 3). The station was located on a uni- 
formly flat surface of crushed rock tailings. No vegetation waf visible 
on or near the dump site. 

Three sets of temperature probes were monitored at this station. One 
set involved a well insulated wooden box filled with Ottawa sand with known 
thermal parameters. Temperature measuremer>ts of the surface and at various 
depths into the sand were recorded during the mission. Two sets of tem- 
perature probes were placed beneath the ground surface. One set was en- 
cased in a "spike" metal tube containing temperature probes at different 
depths. The other set of temperature probes were contained in a plastic 
sewer pipe. These p/'obes were insulated from each other with cotton gauze. 

PRT-5 measurements were made to determine the surface temperatures at 
the station. The field team also recorded the air temperature, near surface 
humidity, wind velocity, cloud cover, net radiation, and radiant flux inci- 
dent to the ground surface for the entire field station. Albedo measurements 
were taken with Exotech (Landsat band) radiometers before and during these 
experiments. Soil samples were also collected at depths corresponding to the 
temperature probes. 
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Two tcts of w«t«r ttmptrature data wart continuously racordad at the 
laach pond northaast of tha Anaconda dump sita (plata 4). Ona proba, 
labalad "A", was placad 10.2 cm from tha bottom of tha pond about } — 
maters away from a drainage pipe where the water flow was nearly constant. 
Tha other proba was placed 12.7 cm from tha pond bottom In a calmer area 
away from any constant flow. Both probes ware mounted In separate wood 
blocks and both were placed about a meter from shore. 

2 . MacArthur Station 

The second recording station was located at the Anaconda Company's 
copper prospect on a wash covered with a sparse cover of low vegetation, 
(Figure 2). Small hills were located on two sides of the site. 

Data similar to that of the Anaconda Dump station was taken at this 
station. A sewer pipe containing insulated temperature probes at various 
depths was used to collect one set of soil temperature profile data. 
Another set of soil temperature profile data was collected by simply 
placing probes into the ground at different depths. The air temperature, 
wind velocity, percent cloud cover, surface temperature, and incident 
short'-waveiength radiant flux were measured every 36 minutes. Albedo 
measurements were also made with the Exotech units. Soil samples were 
collected in order to determine the moisture content profile of the soil. 

B. RESULT.*; 

I . Soil Temperature Measurements 

Temperature profiles for four different times during the 2A hour 
measurement period are found as Figures 3 7* The original data used 
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Plate 4b. LOOKING TOWARDS RECORDING THERMOMETER STATION "B" AT THE 
LEACH POND. 


for thfe graphs are in Appendix 1. The curve.' found for each set cf tem- 
perature probes are more easily described by modelling the temperature 
prof i ie as: 

ilT' ■ AT e where AT' ■ temperature change at damping 

depth 

AT ■ temperature change 
D ■ depth 

DD ■ damping depth 

Thus the depth at which AT' ■ l/eAT is the damping depth. The damping 
depths and temperature variations for the soil at each set temperature 
probe site are tabulated in Table 1. 

The two graphs depicting the temperature profiles at the Macarthur 
site show cooler temperatures for the so!) as depth increases. This is 
a typical summer trend showing the effects of the previous cooling cycle. 

The damping depths are tabulated for each soil site. Both the sur- 
face and a depth Just below it are used to determine the maximum change 
in temperature AT (maximum temperature - minimum temperature) during the 
2k hour period. Results using the surface temperatures are usually better 
to use. However, since the validity of the surface temperature for the 
sand box is doubtful, a near surface depth was also used. AT' ir> then 
c<^<lcu1ated where AT' ■ 1/eAT. By inspecting Figures 3 — 7t one cai deter- 
mine the damping depth by finding the depth where AT for that depth is 
equal to the calculated AT' for the surface or near surface depth. 

The temperature changes during the 2k hour period are shown for each 
soil site at 14. S" and 30" below the surface. The sand box and spike did 
not sense temperatures 30" below the surface. 
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2 . tol*«*M»f r Mddinot 

A ftol'anwur at both ttatlont uMd to tho »hort*««v«Ungth 

radiant flux Incidant to tha ground aurfaca. loth units «nara placad away 
from tha othar Inatrunantt on a ahaat of gray plaatarboard and connactad 
to an ampllfiar and voltmatar to anabla aaay raadout. At tha Anaconda Dump 
aita, tha aol*a-matar waa locatad on flat ground rapraaantat Iva of tha study 
araa. At tha HacArthur alta, a similar aat'up waa placad on a naarby knoll. 

Tha graphad raaults (Piguraa 6 and 9) show that tha radiant flux Inc I* 
dant at HacArthur waa slightly graatar, particularly in tha lata morning. 

t 

This affact could ba tha raault of: calibration problama; a diffaranca in 

tha Illumination or atmoapharic conditions cauaad by haia or clouds; prox> 
imity to raflacting aurfacaa; or problama In lavaling tha units. 

latwaan aunaat and aunriaa (circlad on both graphs) tha maaauramanta 
ara naarly constant, maximum Incidant radlanca occurrad cloaa to 1300 hours 
at both sitas. 

FIguras 8 and 9 show tha plot of sol~a*matar raadings against tlma. 

Tha sol*a*matar maasuras tha radiant flux incidant to tha ground surfaca. 

•2 -t 

Tha original data (Appandlx l)was callbratad to cal cm min amploying 
tha convarsion charts in Appandlx 2. Tha curva for tha uncallbrataH valuas 
is dasignatad by dots, tha callbratad curva by tha "x" symbol. Tha circlad 
points In both flguras danota tha local sunsat and sunrisa. In Figura 8, 
no data was racordad batwaan 2100 and ^<600 hours. 

3. Mat Aadlomatar 

Tha nat radlomatar raadlng raprasants tha diffaranca batwaan tha short 
and long wavalangth radiation hitting a targat and that raflactad back off 
tha targat. Rasults from tha Anaconda dump sita ara glvan In Figura 10. The 
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original data can ba found In Appandin 1. 

Intarpratatlon of tba data Mat compllcatad by tha rounding off of data 
valuat takan Monday Augutt 6. Tha dottad lina through data racordad Augutt 6 
la an approximation aa to haw the curve may have appeared if rounding did not 
occur. Thaaa curvet appear timilar to tha Anaconda dump tita aol'a>metar 
raadingt, hoMOvar, conclutive comparitont ara difficult to tuggett. 

4. Ixotach Landaat land Kadlometera. 

The Cxotech unitt Mare employed to meature the reflectance propertiet 
of the teat aitet. Tmo unitt Mare employed. One unit viewing the target 
with a 1 degree field of v!ew recordt radiant exltance p the aacond unit, 
looking vertically akyward with a 2 pi field of view, racordt total global 
irradiance. The divltlon of the radiant exltance value by global irradiance 
yield the apertured reflectanc/ of the target in the four Landaat wavelength 
Interval t. 

Table 2 littt the albedo value for each tite, calculated by averaging 
the four Landaat reflectance valuet over the wavelength region O.S — 1.1 
micrometert. 


TABLE 2 

Teat Site Albedo 


Site 

Albedo {%) 

Anaconda Dump Site 

22.0 

MacArthur Site 

27.0 

Ottawa Sand 

54.0 


* Some prefer*'bemlapherical conical reflectance". 
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5. Recording Thermometers 

Results from the two recording thermometers at the leech pond ere given 

in Figures 11 end 12. The original date end calibration information are 

located in Appendix 3> Problems were encountered calibrating the two ther> 

«) 

mometers (Appendix 3)* 

Figure 11 gives the thermometer readings every 30 minutes during the 
2k hour data collection period at station A located near the water outlet. 
Figure 12 is a similar graph of station B in a calmer part of the leach pond. 

By comparing these two figures, it is apparent that station A remained 
warmer than station B throughout the 2k hour period. Temperatures at sta- 
tion B varied B.k *C, from 21. S to 29.9 *C, while the variation at station A 
was 7 *C, from 2l».0 to 31.0 *C. The smaller variation at A was anticipated. 
The greater flow of water at A caused greater mixing which would tend to keep 
the water temperature more uniform. However, we anticipated station B to 
have had the highest and lowest recorded temperature, not Just the lowest 
as seen by the data. 

6. Soil Samples 

Soil samples were collected at both stations at depths corresponding to 
those of the temperature probes. Moisture content by percent weight was 
determined for each sample, (Appendix k) . 

Figure 13 shows the moisture contents for the soil surrounding the spike 
and for the soil surrounding the sewer pipe at the Anaconda dump station. 

The points on the curve represent the values determined. 

The moisture content profiles for soils at the MacArthur site are given 
in Figure \k. The dotted line in the sewer pipe curve is based on samples 
whose size was determined to be inadequate to give meaningful results. The 
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DEPTH (INCHES) 

SOIL SURROUNDING SEWER piPB 
rnaep nv»H ^ewp» n.-e ik veht si^au. ^wmk,nts; cess a 


solid linos art the best estimetes of the soil moisture profiles based upon 
the significant date. 

III. THERMAL PARAMETER DETERMINATIONS AND MODELLING 
(AUGUST MISSION) 

The detailed field measurement program carried out in August 1977 at 
Yarington Nevada was designed to supply the surface meteorological and 
thermal data necessary to accurately model the diurnal surface temperature. 
The contemporary principles of one-dimensional heat flow (Carslaw and 
Jaeger, 1959) on the Earth's surface have been converted by a number of 
scientists into mathematical models to resr’ive surface temperature. The 
use and understanding of the modelling results are an extremely valuable 
tool in infrared studies. 

To better understand and analyze the significance of thermal modelling, 
the components of three previously developed models which compute the diurnal 
surface temperature were studied. By comparing the resuits of these models 
using the field data acquired at Yerington, the significance of various 
thermai parameters and the applicability of their rr>sults should be better 
understood. 

The three models compared in this study were: 

1) SURTEMP, a Laplace transform model of the one-dimensional heat con- 
duction equation developed at the Stanford Remote Sensing Lab 
(Lyon, 197^). 

2) WATEMP, a linearized version of the original Laplace transform 
model developed by Watson (1971 and 197^). 

3) CSIROTEMP, a least squares estimation of parameters of surface 
temperature developed by A. Green, (personal comnunication, 1976). 
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Additional models have been developed by Outcalt (1972), Rosema (197^). 
and Kahle (1976), however, they were not available for this study. Previous 
comparisons of the SURTENP and WATEMP models (Marsh, 1975) showad that with 
identical input parameters the model surface temperature results are within 
1 - i|*C. No comparison of the SURTEMP and CSIROTEMP modeis have been previ> 
ously undertaken. 

input parameters for the SURTEMP and WATEMP models can be directly 
extracted or indirectly calculated from the field data. 

A. THERMAL MODEL PARAMETERS 

1) Albedo — the short ‘wavelength (0.5 ‘ 1.1 ym) reflectance of the 
surface as measured by the Exotech radiometers (see Table 2, Section II). 

2) Solar Declination used to calculate solar insolation and deter* 
mined for any particular day from the solar ephemerls, 15.0* for August 8-9, 
1977. 

3) Latitude — the site latitude, necessary to calculate solar Insola- 
tion, 38.9* for the Yeringtor test sites. 

A) Strike and Dip — the orientation of the surface as It Is affected 
by the solar Input. The three Yerlngton sites were essentially level. 

5) Emissivity — the emisslvlty of the site determines the radiation 
temperature observed, most natural surface materials have emissivities in 
the 8 - lA ym region between 0.89 end 0.99* for the modelling work in this 
study the emissivities are approximated at 0.90. 

6) Sky Temperature — the sky back radiation hitting the surface. 

Van Wijk and Shulte (1963) give an empirical formula based upon surface air 
temperature and vapor pressure for the estimation of sky temperature: 





where 


T - T f(P 
sky e w 


T ■ elr temperature 


■ water vapor pressure, and for clear sky conditions 


f(P ) - a. ♦ b.P°‘*5 , a. - 0.678 and b. - 0.0i»l. 
w h h w h h 


Based upon the mean at r temperature and vapor pressure as measured at 
the Yerington Dump site the sky temperature was calculated to be 276. 5*K. 
Measurements made with the net radiometer (Section i 1-6*3) > of net radiation 
(R ) and with the sol-a-meter of short -wave length radiation olong 

with the surface albedo (a), i»«ere employed to determine the sky temperature 


) from the relationship: 


long 


«net 


- 0-.) ♦ 1?'"' ) - "t • 


'short ' '"long ' long' long 
Results from these calculations yield a sky temperature of 272*K. 

Based upon the sensitivity of the SURTEMP enu WATEHP models, an inaccuracy 
of 1 It - 5*K in the value of sky temperature would produce a change in the 
calculated diurnal temperature of 0.S-l.0*C. This preliminarily indicates 
that If a net radiometer is unavailable, the Van Wijk and Shulte (1963) 
empirical relationship should yield reasonably accurate results. 

7) Cloud Cover - Wind Factor — a fractional multiplicative factor 
to compensate for a decrease in solar insolation at the surface due to 
clouds or wind. This factor was estimated to be 0.2 for conditions at the 
Yerington test sites on August 8-9, 1977. 


iaijMbiC 
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B. TEST SITE MODEL RESULTS 


To asteis the eccurecy of thermal inertia modelling reaulta the Stanford 
Remote Sensing Lab has employed (Lyw'> and Harsh, 1976), a standard material 
with known thermal properties in its thermal infrared studies. The standard 
material is 20/30 mesh Ottawa (0.5 -0.8 mm quartz) sand, for which both 
thermal conductivity and heat capacity had been determined by direct labora- 
tory techniques at varying moisture contents (Hoench, 1969). 

The sand was placed in a well insulated (plate I) Douglas fir wood box 

at the Anaconda Dump site. The moisture content of the sand was approxi- 

33 - 2-1 

mated at 0.02 cm /cm which wou.J give a thermal Inertia of 0.03 cal cm *C 

sec The SURTEMP and WATEMP models were then run to determine the accuracy 

of the input parameters. By comparing the observed surface temperatures of 

the sand, (recorded by the PRT-5) with the mode* predicted temperatures, for 

20 equally spaced time increments in a 24 hour period, the error limits r'f 

the input parameters and models can be analyzed. Table 3 gives the input 

parameters for the Ottawa sand, and Table 4, the standard error of the fit 

(SE) of the diurnal surface temperature ^or a range of thermal inertia 

values. The standard error of the fit is calculated from the formula: 


SE - V*/"-’ ’7^)' 

1 me 


where n 


T 

m 

T 

c 


the number of sample points (20) during a 24 hour 
period 

measured surface temperature of the Ottawa sand 

model calculated surface temperature of the Ottawa sand. 
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MODEL PARAMETER VALUES FOR YERIN6T0N SITES. AUGUST 8,9 197 7 


Anaconda Dump Sitf McArthur Sitt 

Cruahadjaachad Granodlorlta solli 

quartz monzonita 


Albado 

22.0 

27.0 


Emissivlty 

0.90 

0.90 

Cloud Covar 

0.2 

0.2 


Latitude 

38.9 

38.9 


Declination 

15.0 

15.0 


Dip 

0.0 

0.0 


Strike . 

0.0 

0.0 


Sky"Tamperature" 272K 

272K 


Moisture Content 

Solke Sawer pipe 

-1" ’TT* iTo 

-10" 4.5 3.4 

-16" 6.0 3.5 

Spike 

1.0 

Sewer pipe 
1.3 
2.5 
3.0 
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TABLE 3 

Input Parameters for the Ottawa Sand Study 


Albedo : 

0.S4 

Solar 

OecI inaticMi 

15.0 

Emissivlty : 

0.90 


Lat 1 tude 

38.9 

Cloud Cover : 

0.20 


T . 
sky 

272*K 

Dip : 

0.0 

TABLE A 

Thermal Model Accuracy 

Strike 

0.0 

T. 1 . 


Moisture Content 

Standard Error of Fit 

-2.-1 -iv 

cm C sec ) 


(cm^/cm^) 

SURTEMP 

WATEMP 

0.017 


0.0 

6.66 

11.15 

0.030 


0.02 

5. AS 

9.2A 

0.037 


O.OA 

5.65 

9.25 

0.0A3 


0.09 

5.96 

9 .A 9 


Both models produce the minimum error between observed end model pre** 

•2 “I “i 

dieted surface temperature with a thermal inertia of 0.03 cal cm *C sec 
The SURTEMP model best fit produces an error of about 5.5*C, the WATEMP model 
9.2li*C. These results indicate that both models appear to correctly approxi* 

t 

mate the thermal inertia of the Ottawa sand under the meteorological condi- 
tions present August 8-9, 1977. However, their ability to reproduce the 
surface temperature is at best within 5 - I0*C. 


Tht CSIROTEHP fflodcl rtquircs Input of the »ltc latitude, albedo, tolar 
declination, and observed temperature through the diurnal cycle. Estimates 
of the mean surface temperature, thermal Inertia, and back radiation are also 
Input. The model then revest Imates these terms to arrive at a least squares 
estimate of these values based upon the observed surface temperatures. In 
an alternate mode the mean surface temperature and back radiation arc re* 
estimated keeping the thermal Inertia at the original input value. Model 
results are given in Table 5 for the expected range of thermal inertia 
values kept constant. 


TABLE 5 

CSIROTEMP Model Results 


T.l. 

(cal cm ^ *C 'sec 

0.020 

0.030 

O.OkO 


Standard Error of Fit 
CSIROTEMP 

3.1 

1.8 

1.3 


CSIROTEMP model results allowing the thermal inertia to be re~estimated 
to arrive at a best fit calculated the thermal inertia of the Ottawa sand to 
be 0.037 cal cm*^ *c"' sec"*. This is off by 0.007 from the SURTEMP and 
WATEMP results, however. It Is important to note the extremely small differ- 
ence in the standard error between a thermal inertia of 0.030 and 0.037 in 
all three models. It appears reasonable to conclude from these results that 
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FIGURE 15a Yerington Dump Site .Sandbox Experiwent, WOE 1 —all estiwated. 
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FIGURE 15b Yerington Dunp Site, Sandbox Experiment— CSIR0TEW> Ave. tei^ and Thermal 
Inertia estimated. 






parameter (G) to be estimated. 






all three models were able to calculate the thermal Inertia of the Ottawa 
sand within approximately 20% based upon the input parameters and their 
associated error limits. 

The next step in the analysis employed the SURTEMP model to estimate 
the thermal Inertia of the Anaconda Dump and MacArthur test sites. The 
input parameters used in the Ottawa sand run were kept constent, changing 
only the albedo of the site as determined from the Exotech measurements. 
The stanriard error of the fit between the observed and model calculated 
diurnal temperatures for the anticipated range of thermal Inertia values 
were then calculated (Table 6). 

TABLE 6 

Yerington Test Site SURTE.MP Model Results 


T .1 . 

Standard 

Error of Fit 

(cal cm ^ *C ' sec 

Anaconda Dump 

MacArthur 

0.020 

^ mm m 

6.57 

0.030 

12.36 

5.76* 

0.035 

11.92 

6.2A 

O.OkO 

11.69 

6.87 

0.0A5 

11.58* 

l.kl 

0.050 

11.59 

— 

0.060 

11.62 

— 

* 

■ best fit 
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These results indicate the thermal inertia of the test sites can be 
estimated within 10 to 15 percent at 0.030 for the MacArthur site and 0.045 
for the Anaconda Dump site. The error at the Anaconda Dump site in the fit 
of observed to calculated diurnal temperature is nearly double that of the 
HacArthur and Ottawa sand sites. This fact is difficult to explain unless 
the albedo or local meteorology was more variable than evidenced in the 
measurements. 


i 


C. APPROXIMATION TO A RADIANT HEATING METHOD FOR OFTERMINING THERMAL 
INERTIA 


Schultz ( 1968 ) developed a novel method for nondestructively determining 
the thermal inertia of solids near ambient temperature. The method Involved 
heating, with radiant energy, a small area on the surface of a *'semi -inf in i te" 
solid for a short period of time. The characteristically shaped temperature 
rise recorded with an IR radiometer is compared with that of a reference 
standard yielding the thermal inertia of the sample from the relationship: 


As a means of approximating Schultz's method the Ottawa sand was 
employed as the reference standard and the natural solar heating was em- 
ployed, over a defined period of time, as the source of radiant energy. 

An early morning time period (0710-0825) was chosen to run the experiment 
due to the more stable meteorological conditions and faster heating in this 

time frame. Results from the calculations yield a thermal inertia of the 

-2 -1 -i ~2 - 

MacArthur site of 0.029 cal cm C sec and a value of 0.031 cal cm C 

.4 

sec for the Anaconda dump site. 


4 

i 
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The MacArthur thermal inertia result is quite close to the value deter* 
mined by the modelling, however, the Anaconda Dump site is off by nearly 35^* 
From these results It Is Impossible to claim success for this approximation 
to Schultz's method. The numerious influencing meteorological variables 
which existed during our heat inn history obviously detracted from the use* 
fulness of our technique. Under Schultz's ideal conditions there are no 
unknown variables and the heating rate 's kept constant. The development 
of this method Into a field instrument with Its own heat source might be a 
very useful next step in thermal lner".Ia studies. 

IV. SECONDARY FIELD STUDY - YERINGTON NEVADA 
(DECEMBER 8-^, li)78) 

A second field mission was carried out coincident with the December 8th 
and 9th (1978) HCMM satellite overpass of the Yerington area. Due to weather 
conditions and availability of personnel, it was impossible to complete as 
detailed a measurement program as accomplished during the August 1977 mission. 
Data recorded during this mission Included ground surface temperatures (PRT-5), 
subsurface soli temperatures at i, 2, and 8 inches (thermistor probes), air 
temperature, windspeed, and cloud cover conditions. 

The sites chosen for this study were designed to represent both hydro* 
thermally altered (MacArthur) and unaltered (Mason Butte) areas, as well as 
one site in Alluvium. 

A. DESCRIPTION OF TEST SITES 

Mason Butte (I): a small butte (1.4 x 2.8 km) (plate &) north of the 

Yerington town site. Test site 1 Is an area of unaltered Jurassic granodiorite 
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pebble~cobble s}il surrounded by greesewood end shedescale phreatopHyte 
vegetation. 

Hason Butte (2): (plate 7) a small Jagged outcrop (2 x 5 m) of Jurassic 

unaltered granodiorlte. 

Alluvium (3): (plate 8) a medium to fine grain brown to tan soil 

derived from the Tertiary ignimbrite sequence west of the test site. The 
vegetation at the site is predominantly greasewood and sagebrush. 

HacAr':hur {k): (plate 6) undisturbed soil at the HacArthur prospect 

an area of hydrothermal ly altered granodiorite and quartz monzonite. 

HacArthur (5): (plate 8) an unvegetated flat mound created by trenching 

of the area of hydrothermal I y altered granodiorite, quartz monzonite porphyry, 
and andesite cobbles. 

B. THERMAL DATA 

Data collected for the five test sites are given In Table 7, all temper- 
ature values are given in degrees C. A plot of the diurnal temperature pro- 
file is given in Figure 12. 

In an attempt to establish a calibration site large enough to be appli- 
cable to the HCMM system, data was collected from the large (0.5x 1.5 km) 
cooling ponds for the Fort Churchill Power Station. The power station is 
located approximately 15 km north of the Yerington test sites in Mason Valley. 
Plant overflow and inflow water temperatures were supplied for the approximate 
times of sateliite overpass, and are given in Table 8. 
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TABLE 7 

Thermal Data Sheets — Yerington — December 8-9, 1978 


Date 

Time 

PRT-5 

Therm 

i" 

Probe i" 

Probe 2*' 

Probe 8" 

Air 

Temp 

Wind 

Cloud 

Mason Butte Soil (1) 







12/8/78 

0812 

•12.0 


<0 

<0 

<0 

-4 

0-2 mph 
High 10% 


09^5 

- 4.0 

-5.5 

<0 

<0 

<0 

-6 

0-T 

20% 


1115 

♦ 3.5 

-0.5 

<0 

<0 

<0 

-1 

0-2 

30% 


1328 

3.5 

4.2.4 

<0 

<0 

<0 

+0.5 

0-2 

50-60% 


1403 

♦ 5.0 


2.0 

<0 

<0 

-0.5 

0-2 

80-90% 


1535 

0.0 


<0 

<0 

<0 

-0.5 

0-2 

90-100% 

Sunset 

1619 

- 1.5 


<0 

<0 

<0 

-3.5 

0-2 

80-90% 


1725 

- 3.0 


<0 

<0 

<0 

-4.5 

0-2 

40-80% 


1735 

- 3.5 


<0 

<0 

<0 

-4.5 

0-2 

40-80% 

12/9/78 

0030 

-11.5 


<0 

<0 

+0.2 


0 

0 


0155 

-10.: 


<0 

<0 

+0.8 


0-2 

0-10 


0209 

-10.5 


<0 

<0 

+0.2 


0-2 

0-10 


0915 

- 3.0 


<0 

<0 

<0 


0-2 

30-40% 


1110 

♦ 7.5 


+2.5 

<0 

<0 


0-2 

20% 


1201 

4 7.5 


+4.2 

<0 

<0 


0-2 

20-40% 


1335 

♦ 9.0 


+6.1 

<0 

<0 


0-2 

50-80% 


Ha son Butte Outcrop (2) 

12/8/78 0812 -11 
09A5 - 5 
1115 ♦ 0.5 
1328 ♦ 5.5 
1A03 * 2.0 
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TABt»I 7 (contlnu*d) 


r 


Date 

Time PRT-5 Probe i" Probe 2" Probe 8" 

Air 

Tamp 

Wind 

Cloud 



Hason Butte Outcrop (2) (cont'd) 


12/8/78 

1535 

0.0 




1619 

- 2.0 




1725 

- 3.5 




1735 

" 3.5 



12/9/79 

0030 

-10.0 




0155 

1 

O 




0915 

- 5.0 




1110 

4 2.0 




1201 

♦ A. 5 




1335 

♦ 8.0 



Alluvium 

Ha son 

Pass (3) 



12/8/78 

0836 

- 5.5 -8 

<0 

<0 


1005 

♦ 2.0 -2 

<0 

<0 


1100 

4 8.0 0 

<0 

<0 


1310 

4 6.5 41. A 

41. A 

<0 


Ili20 

4 5.5 

41.2 

<0 


1517 

4 1.0 

40.1 

<0 


1635 

- 0.5 

0.0 

<0 


1710 

- 2.5 

<0 

<0 


1750 

- 5.5 

<0 

<0 

12/9/79 

OOkS 

-10.0 

<0 

<0 


0137 

- 9.0 

<0 

<0 


<0 

<0 

<0 

<0 

<0 

<0 

<0 


<0 

<0 


<0 


- 5 

- k 
0 

♦ 0.5 
0 

♦ 2.0 


- 3.0 


- 2.0 


" 2.5 


- 13.0 

-12.5 


0-2 mph 
10 % -- 
0-1 

20 % 

0-2 

20 % 

0-2 

AO-50% 

0-2 

80 - 90 % 

0-1 

90-95% 

0-2 

90 - 100 % 

0-1 

AO-80% 

0-2 

50% 

0 

0 

0-2 

0 - 10 % 


<0 






50 


TABLE 7 (eontinu«d) 


0at« 

Tl 

PRT-5 

Tharm 

i" 

Proba i" 

Proba 2" 

Proba 8'* 

Air 

Tamp 

Wind 

Cloud 

Alluvium Haton 

Pa»» (3) (cont'd) 





12/9/79 

022A 

-10.0 


<0* 

<0* 

<0* 

-13.5 

0-2 mph 
0-10!fe 


0935 

♦ 5.5 


<0 

<0 

<0 

0.0 

0-2 

yo% 


1055 

♦ 9.5 


♦ 0.2 

<0 

<0 

♦ 1.5 

0-2 

30!L 


1222 

♦Ili.O 


♦ 2.8 

<0 

<0 

♦ 5.5 

0-2 

25% 


1313 

♦16.5 


♦ 3.6 

<0 

<0 

♦ 8.0 

0-2 

60-50% 

MacArthur Soil 

(6): 







12/8/79 

0851 

- 6.0 

-5 

<0 

<0 

<0 

- 3.0 

O 

1 

o 


1021 

♦ 5.0 

♦1 

<0 

<0 

<0 

- 1.0 

0-1 

20% 


I0l»7 

♦ 8.0 

-0.5 

0.0 

<0 

<0 

- 3.0 

0-1 

20% 


1300 

♦ 9.0 


♦ 3.5 

<0 

<0 

♦ 0.5 

O 

1 

o 


1A30 

♦ 5.0 


♦ 1.5 

<0 

<0 

- 1.5 

0-2 

80-90% 


1508 

♦ 0.5 


♦ 0.2 

<0 

<0 

- 2.5 

0-2 

90% 


16A5 

> 1.5 


<0 

<0 

<0 

- 1.0 

0-2 

90% 


1703 

- 2.C 


<0 

<0 

♦0.1 

- 5.0 

0-2 

90% 


1800 

- 6.0 


<0 

<0 

♦0.1 


0 

1 

o 

M/S/7S 

0058 

- 8.0 


<0 

<0 

♦0.3 


0 

c 


0I2A 

- 8.0 


<0 

<0 

♦0.2 


0-2 

0-10% 


0238 

- 9.0 


<0 

<0 

♦0.2 


0-2 

0-10% 


0950 

♦ 2.5 


0.0 

<0 

<0 


0-2 

60% 


lOAl 

♦ 5.0 


♦ 2.6 

<0 

<0 


0-2 

30% 


1235 

♦ 8.5 


♦ 6.3 

<0 

<0 


0-2 

25% 


1305 

♦ 7.5 


♦ 6.7 

<0 

<0 


0-2 

60-50% 
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TABLE 7 (continuid) 


Dat« 


Tim# PRT-5 Prob# i" Prob# 2" Prob# 8" 


Wind 

Cloud 


MacArthur Mound (5) : 

12/8/79 0851 - 7.0 
1021 • 1.0 
10A7 * 2.0 

1300 4 3.0 

1A30 4 1.5 
1508 4 0.5 
16A5 - 2.0 
1703 - 2.0 
1800 * it.O 

12/9/79 0058 - 9.5 
012A • 9.5 
0238 -10.0 
0950 0.0 
lOAl 4 5.0 

1235 4 6.5 
1305 ♦ 5.5 


VeRiN&TOil RCHU 
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OF POOR QUALITV 



9L/klei €i/«/ei 

0090 OOhO OOOO 0001 0001 oott 
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TABLE 8 

Fort ChurchMt Power Station Cooling Pond Temperatures 


TIME 

OUTFLOW (*C) 

INFLOW (*C) 

MEAN (*C) 



0200 

23.3 

12.8 

18.0 

1300 

31.7 

13.9 

22.8 


C. MODEL RESULTS 

Unfavorable weather conditions and a severe lack of time made it 
Impossible to transport the Ottawa sand standard to the field during this 
mission. Therefore, the meteorological input variables could only be 
estimated. Cloud cover was visually estimated during the measurement 
program and day and night sky temperatures were evaluated employing the 
Van Wijk and Shulte equation. The albedo of aach site was determined from 
previous measurements of the areas with the Exotech radiometers. 

The SURTEMP model was run to estimate the thermal inertia of the 
Mason Butte and MacArthur test sites employing the input variables given 
in Table 9. 

The standard error of the fit between the observed and model calcu* 

-2 

lated diurnal temperatures for the anticipated range (0.020-0.060 cal cm 

• 1 * 1/2 * 

*C sec ) of thermal inertia values was then calculated. The results 

indicated the thermal inertia of the test sites could be est^imated with 

with valines of 0.035 for the MacArthur site and 0.050 for the Mason 

Butte site. The standard error of the fit for the MacArthur site is 
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rees and is 3.96 degrees 

for the Mason 

Butte site. 


TABLE 9 


Input Variables for 

the December 

1978 Mission 


Mason Butte 

MacArthur 

Albedo 

0.22 

0.27 

Emissivlty 

0.90 

0.90 

Cloud Cover 

0.50 

0.50 

Latitude 

38.90 

38.90 

Solar Declination 

22.70 

22.70 

Dip 

0.0 

0.0 

Strike 

0.0 

0.0 

Sky Temperature 

220. OK 

220. OK 


A value of 0^0^ for the granodiorite at Mason Butte Is considered 
reasonable In light of previously published (Lyon, I97it) values for similar 
rock compositions. The value of 0.035 for the altered granodiorite and 
quartz monzonite soil is within 15% of the value determined during the 
August 1977 field mission. This is considered to be within experimental 
and model limitations due to the quite variable meteorological conditions 
and the associated uncertainty of relat fg these conditions to model input 
variables during this mission. 








MASON BUTTE SOIL, TERINGTON 



FIGURE 15d. Yerington Area, Mason Butte Soil, Mode 2— CSIROTEMP Ave.Temp 
and Thermal Inertia to be estlemted. 





MASON OUTCROP (2) , YERINGTON 

-KOOE 1 VBflR.P RND TSKT (Gl RLL TO BE E5T1MHTED 
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(3 030) 3UnibU3dN3i 








FIGURE 15e Yerington Area .Mason Butte. Granodiorite Outcrop. CSIROTEIC 

Mode 1— Ave Temp. Thermal Inertia, and Atmospheric Par^ter (G) 
to be estimated. 
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V. P-3 (MMS), U-2, AND HCMM DATA INTERPRETATIONS 


The chart below describes the data available for the study: 

August 8-9. 1977 May 30. 1978 December 8-9. 1978 


Field Measurement 
P-3 (MMS) 

U-2 

HCMM 


X 

X 

X X 

X 


X 


X 


The field measurements have been described elsewhere. Following will 
be an interpretation of the P-3 (MMS), U-2, and HCMM data. 


P-3 (MMS) DATA 

The P-3 (MMS) flight over Yerington was made on August 8th and 9th 
to correlate overflight data with three local ground calibration sites. 

The data obtained by the P-3 (MMS) corresponds to 11 channels ranging from 
O.A 19 to 11. 6A pm (channel 11: 7.886 - 11. 6A ym) . The digital data re- 

corded in magnetic tapes was analyzed to estimate the earth surface tempera- 
tures. 


Method to Calculate Video Pixels Temperatures 

The equation to convert the (pixel) counts to temperature is done using 
the formula provided by NASA/JSC 


^BB(Hi) ~ ^BB(Lo) 
HKW(Hi) - 


X 


n 


'"'“(Lo) 


f 
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h 

« 


T 

n 

■ 

Temperature in *C o* pixel n 


^BB(Hi) 

m 

Temperature 'lack Body - High 


"^BB(Lo) 

m 

Temperature Black Body — Low 


HKW(hi) 

m 

Binary Count in Housekeeping Word for 

Black Body — High 

«"^Lo) 

m 

Binary Count in Housekeeping Iford for 

Black Body ~ Low 

X 

n 

m 

Binary Count of pixel n. 



The digital data was e>;tracted and analyzed with the following 
resul ts : 

1) For the day flight (tape 000162, Mission 366) the low and high 
temperatures for the black bodies were reversed, which means 
that the data were incorrect and probably useless. The values 
encountered were: 

Housekeeping Word 923 * 21 (Low Temp 6B1) 

Housekeeping Word 926 ■ 18 (High Temp BB2). 

2) For the night flight (tape OOOiil , Mission 366) the range of 
temperatures for channel 11 were too widely spread to use in 
the calculation of pixel temperatures: 

Housekeeping Word 923 " 26 (Low Temp BBl) 

Housekeeping Word 926 ■ 251 (High Temp BB2). 

Deu to the problems described above, it was not possible to obtain meaning- 
ful temperatures from the P-3 (MMS) flight to correlate with the ground 
measurements. 

U-2 (HCMR) DATA 

Two U-2 aircraft missions were conducted over the test site to corre- 
late overflight data with ground calibration sites. The first mission was 
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flown In August 8-9, 1977i coincident with the P-3 (MMS) flight. The data 
obtained by the Heat Capacity Mapping Radiometer corresponds to two channels 
(v’sible: 0.5 - 0.7 Um, thermal Infrared: 10.5 - 12.5 pm) with a nominal 
thermal resolution of 0.2*C and spatial resolution of 57><57 m. 

1 . August 6-9, 1977 Flight 

Two tapes were received for the August 8-9, 1977 flight: U2-HCMR- 

ViCAR and U2-HCMR-cal ibrated. The analysis was performed on the U2-HCMR 
calibrated data. Figure 13 is a "panoramic corrected" OOTPRINT showing 
part of the Weed Heights near Yerington, Nevada. The OOTPRINT was obtained 
using the raw numeric data in the thermal infrared channel from the \j Z 
HCHR calibrated tape (U-2 Flight 77.130, HCM Flight l8). The area inside 
the square has been enlarged and is shown in Figure 14. The raw numeric 
data (channe' 2) corresponding to the area inside the square is shown in 
Table 10. The numbers in the figures represent the coordinates of the U2- 
HCMR calibrated tape. The enlarged area in Figure lit is part of a water 
body (tailing pond) in the Yerington mine dump area, for which there is 
ground truth data. Note that the image (Figures 13 and 1li) has been 
"reflipped horizontally", i.e., rotated from how it was originally taken, 
and now west is to the left, and east to the right. 

In order to analyze the numeric data corresponding to the radiance 
and temperature for channel 2 and the reflectance for channel 1 (visible), 
data were extracted from the U2-HCMR calibrated tape corresponding to the 
area shown in Figure 1A and compared with the IBM output provided by 
NASA/GSFC. Table 10 lists the raw numerics and calibration data from the 
U2-HCHR calibrated tape, and Table 11 shows the calibration data from the 
IBM output. Tables 10 and 11 indicate that there is no correspondence 
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bttwten tht»e radiance and fmptfture sets , which m«/ns thit some problems 
(not yet clarified) exist In the data. The table, he ;ver, does show that 
there is agreement between the IBM output and the U2-HCMR calibrated tape 
for the raflactanca data In the visible . 

The chart below compares the temperatures between the U2-HCMR cali- 
brated tape and the ground measured temperatures, in degrees Kelvin: 


STATION 

U-2 TEMPERATURE (*K) 

GROUND TEMPERATURE (*K) 

A (water body) 

NASA/GSFC 

303.28 

IBM Output 
304.97 

295.65 (Hr: l4:4l) 

B (water body) 

312.65 

317.01 

296.05 (Hr: 14:45) 

AT 

9.37 

12.04 

0.40 


The M for the U-2 data (over water) is too large considering that stations 
A and B are one pixel apart (see Figure 14). For the Dump Area station the 
U2-HCMR calibrated tape temperature is 320.86 (NASA/GSFC) and the ground 
measurement 306.2 (AT ■ 14.66). Figure 15 shows a DOTPRINT with the loca- 
tion of the dump area, and Table 12 gives the data extracted from the U2- 
HCMR calibrated tape. 

2. May 31. 1978 

The Hay 31, 1978 mission was conducted by NASA/ARC without informing 
us until 6 weeks after. Needless to say we have no ground measurements. 

No analysis was performed on the May 31, 1978 mission over the test 
site, due to the lack of correct calibrations for the data and ambigu ity 
of the data. 
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1 



Figure Enlarged "panoranlc corrected" Dotprlnt for part 
of the ‘tailing pond. The raw numeric data for the infrared 
channel from Table lO(U2-HCMR) was used to produce the Dot- 
print • ... — - 


\ 





b3 
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Figure 15. ENLARGED "PANORAMICALLY-CORRECTED DOTPRINl 
FOR THE DUMT AREA (1)1) 


TABLE 10. Raw Numaric and Calibration Data Extracted from tha U2>HCMf. 

Calibrated Tape 


HCH HEAT CAR. MAPPER 
U-2 FLIGHT 77.130 
HCH FLIGHT I8 
YERINGT0.I. NEVADA 

HCMR FLIGHT INSTRUMENT MASTER CALIBRATION TABLES 

AREA IS 576 TO 58 1 SCANLINES DOWN FROM THE START-0 
1A5 TO IA9 PIXELS IN FROM LEFT EDGE-0 


INFRARED BAND 


RAW DATA MATRIX 

(WINDOW) 




COLUMN: 

(165) 




(169) 

ROW: (576) 

500 

512 

531 

660 

605 


501 

538 

665 

668 

623 


506 

535 

666 

668 

536 


533 

666 

667 

682 

660 


502 

656 

670 

591 

625 

(581) 

651 

667 

588 

639 

628 


CALI 

1 BRAT ION TABLE 

(EXTRACT) 




RAW 

BYTE# 

OCTAL 

RADIANCE 

OCTAL 

TEMPERATURE 

665 

1777:2000 

7637366552 

.00192150 

10306563065 

302.388 

666 

1781:1786 

7637606102 

.00192^10 

10306563653 

302.668 

651 

1081 ... 


.00192110 


302.766 

656 

1813 ... 


.00193590 


302.926 

660 

1837 


.00196550 


303.283 

66L 

1853 


.00195190 


303.521 

667 

1865 


.00195670 


303.699 

668 

1869 


.00195830 


303.758 

670 

1877 


.00196150 


303.877 

682 

1925 


.00198070 


306.586 

500 

1997 


.00200950 


305.662 

501 

2001 


.00201110 


305.700 

502 

2005 


.00201270 


305.759 

506 

2013 


.00201590 


305.876 

512 

2065 


.00202870 


306.362 

531 

2121 


.00205910 


307.663 

533 

2129 


.00206230 


307.558 

535 

2137 


.00206550 


307.673 

536 

2161 


.00206710 


307.731 

538 

2169 


.00207030 


307.866 

588 

2369 


.00215030 


310.696 

591 

2361 


.00215510 


310.863 

605 

2617 


.00217750 


311.669 

623 

2689 


.00220630 


312.653 

625 

2697 


.00220950 


312.765 
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TABLE 10. (continued) 


CALIBRATION TABLE (EXTRACT) 


RAW 

BYTE# 

RADIANCE 

TEMPERATURE 

628 

2509 

.002211(30 

312.932 

639 

2553 

.00223189 

313.5<»1 

6A0 

2557 

.002233^9 

313.596 


VISIBLE BAND 

ftAW 5AtA Matrix (window) 

COLUMN: (1A5) 0^9) 

ROW: (576) 2k2 2$U 298 250 307 

2<«6 303 250 2A6 307 

263 311 250 2A2 268 

311 263 2A6 250 316 

285 250 2A2 290 368 

(581) 255 2A6 29^* 350 398 


CALIBRATION TABLE (EXTRACT) 


RAW 

BYTE# 

REFLECTANCE 

2A2 

965 

.231 

2A6 

981 

.235 

250 

997 

.239 

255 

1017 

.2l(i( 

263 

10i(9 

.252 

268 

1069 

.257 

285 

* 1137 

.21k 

290 

1157 

.279 

29^ 

1173 

.283 

303 

1209 

.292 

307 

1225 

.296 

311 

12M 

.300 

316 

1261 

.305 

350 

1397 

.339 

368 

1A69 

.357 

398 

1589 

.387 


T 
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TABLE 11. Calibration Data from IBH Output (NASA/CSFC) 


HCMR FlICHT MASTER CALIBRATION TABLES 
NAME- §76 jI'089 
LENGTH- 102A ELEMENTS 



TABLE GENERATION PARAMETERS 



INFRARED 

VISIBLE 



BASE INCREMENT 

BASE 

INCREMENT 

LOCATION 

VALUE 

LOCATION VALUE 


2 

1.029itE-03 2.1000E-06 

11 0.0 

1.0000E-03 


INFRARED 

JWND 


VISIBLE BAND 

INDEX 

RADIANCE 

TEMPERATURE 

INDEX 

REFLECTANCE 

kkS 

.0019597 

303.8103 

242 

0.231 

kke 

.0019618 

303.8879 

246 

0.235 

A51 

.0019723 

304.2761 

250 

0.239 

i»54 

.0019786 

304.5083 

255 

0.244 

460 

.0019912 

304.9717 

263 

0.252 

464 

.0019996 

305.2798 

268 

0.257 

467 

.0020059 

305.5105 

265 

0.274 

468 

.00200C0 

305.5874 

290 

0.279 

470 

.0020122 

305.7407 

294 

0.283 

482 

.0020374 

306.6577 

303 

0.292 

500 

.0020752 

308.0220 

307 

0.296 

501 

.0020773 

308.0974 

311 

0.300 

502 

.0020794 

308.1726 

316 

0.305 

504 

.0020836 

308. 3232 

350 

0.339 

512 

.0021004 

308.9241 

368 

0.357 

531 

.0021403 

310.3408 

398 

0.387 

533 

.0021445 

310.4890 



535 

.0021487 

310.6372 



536 

.0021508 

310.7112 



538 

.0021550 

310.8591 



588 

.0022600 

314.5090 



591 

.0022663 

314.7251 



605 

.0022957 

315.7295 



623 

.0023335 

317.0110 



625 

.0023377 

317.1526 



628 

.0023440 

317.3650 



639 

.0023671 

318.1411 



640 

.0023692 

318.2114 




i- 

k 
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TABLE 12. U2>HCHR Calibrated Tape, Second Version, Fit 18 (Day) 5 by S 
Window, Calibration Site it Approximately at the Center of the 
Window or Right or Left One Pixel. (Possibly Evan in Any of 
Middle 3 Rows, Any of 5 Pixels) 

Note, All Row Coordinates are From Top'Line at Zero. All Pixel 
Coordinates are From Lcft'PIxal at One, as Per Documentation 
of Tape Format. (Indicated by Upper Left Corner«ULC) 


INFRARED CHANNEL 

FORMAT; TOP-RAW NEXT-INDEX NEXT-RADIANCE BOTTOM-TEMPERATURE 



(I)!'*? 

148 

149 

150 

151 

(0) 

ULC 





589 

705 

747 

730 

721 

692 raw • 


2820 

2988 

2920 

2884 

2768 index 


.00250570 

.00259390 

.00255820 

.00253929 

.00247840 radiance 


322.71655 

325.56055 

324.41553 

323.80615 

321.82593 

590 

707 

647 

661 

645 

747 


2828 

2588 

2644 

2580 

2988 


.00250989 

.00238390 

.00241330 

.00237970 

.00259390 


322.85327 

318.70288 

319.68140 

318.56274 

325.56055 

591 

705 

685 

678 

659 

667 


2820 

2740 

2712 

2636 

2688 


.00250570 

.00246370 

.00244900 

.00240910 

.00242590 


322.71655 

321.34424 

310.86108 

319.54199 

320.09863 

592 

666 

705 

717 

692 

685 


266i« 

2820 

2868 

2768 

2740 


.00242380 

.00250570 

.00253090 

.00247840 

.00246370 


320.02930 

322.71655 

323.53442 

321.82593 

321.34424 

593 

690 

726 

695 

731 

717 


2760 

2904 

2780 

2924 

2868 


.00247420 

.00254980 

.00246470 

.00256030 

.00253090 


321.68848 

324.14502 

322.03198 

324.03198 

323.53442 




VISIBLE CHANNEL 



FORMAT: TOP-RAW NEXT- INDEX BOTTOM-REFLECTANCE 


(1)147 

148 

149 

150 

151 

(0) 

ULC 





589 

385 

363 

359 

359 

372 


1540 

1452 

1436 

1436 

1488 


.374 

.352 

.348 

.348 

.361 

590 

372 

389 

398 

385 

372 


1488 

1556 

1592 

1540 

1488 


.361 

.378 

.387 

.374 

.361 

591 

398 

394 

394 

389 

402 


1592 

1576 

1576 

1556 

1608 


.387 

.383 

.383 

.378 

.391 





592 389 

1556 
.378 


381 

152<« 

.370 


TABLE 12 (continued) 


593 385 372 

1SA0 1688 

.376 .361 


363 

372 

389 

1632 

1686 

1556 

.352 

.361 

.378 

381 

368 

363 

1526 

1672 

1632 

.370 

.357 

.352 


Oh<t:N ■ o*- •'» 
OF T'OOH QUALITY 


69 


H CMM DATA 

The mein reason to use HCMM data In the present study is to try to 
discriminate different rock types using thermal parameters, parcicularly 
those related to rock density. The HCMM sensor is a two>channe1 scanning 
radiometer providing measurements of reflected solar (0. >5 to 1.1 ym) and 
emitted thermal energy (10.5 to 12.5 ym) . There are satellite night/day 
coverage patterns at least once every M days at approximately 12*hour 
intervals. At northern hemisphere mid -latitudes, the crossing times are 
about 1:30 p.m. and 2:30 a.m. From the nominal orbit altitude of 620 km, 
the spatial resolution of the infrared channel is approximately 600 by 
600 meters at nadir, and the resolution in the reflectance channel is 500 
by 500 meters. These values are masked ly resampling the data in the 
data processing, which generates registered data at a 1<3l.5 meter pixei 
size . 

The data obtained by the HCMM are digitized to 8 bits units of energy 
(255 levels). In this range, the byte (pixel) counts can be converted to 
temperature using a formula derived from the Planck function and from 
calibration procedures representing radiometer performance. The pixel 
counts are converted to temperatures using the formula of Bohse et al., 
1979 established according to a performance evaluation of the HCMM: 

T(l) - Kjj/An^K^/l - l] 

I « image pixel data value (0-255) 

T ■ temperature in degrees Kelvin 

- 11|k21.587 
Kj • 1251.15191 

Kj • -118.21378 
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For • given I value (pixel count) we believe we can calculate the tempera* 
ture of a surface element in *K (degrees Kelvin) for the range: 

1-0 T - 260*K t • -13*C 

1- 255 T - 340*K t - 467*C 


Interpretation of HCMM*Satel 1 1 te Registered Sun^ e Temperatures for 
30th MI.V. 1978 Night IR (scene A*A003^*102t0) . 


The HCMM night* I R data for May 30, 1978 was analyzed for the study 
area. The analysis was conducted using (1) the 3*‘Color Grinnell TV*di splay 
which allows viewing of selected ranges of digital data to construct false* 
color composites of the IR channel, and (2) simulated grey*level pictures 
(DOTPR.'NT) on a Printronix matrix printer. 

Figure 16 is a DOTPRINT at scale 1:250,000 where the cold areas are 
represented in dark and warm areas in white with grey-leveis In between. 

The test*site corresponding to the Yerington mine is clearly identified 
In the OOTPRINT with the symbol Y O 


Figure 17 is a densi ty*sl iced DOT 
increment of 0.38^C between each class 
sented. All the rest (cool areas) is < 
ture for the terrain in the Yerington i 
Range In *K 

281.08 • 280.70 
279.56 - 279.18 
278.80 - 278.42 


lINT at 1:250,000 scale with an 
There are five classes repre* 
lite. The differences in tempera* 
ne are from warmer to cooler: 
Terrain 

Tailing Ponds — water 
Rock Dumps 
Waste rock 


Figure 18 is a density-sliced DOTPRINT at 1:250,000 scale. Class A 
represents temperatures from 287.5 *K to 281.8 *K and correspond to the 
warmer areas (water). Class B represent temperatures from 281.4 *K to 


1.-00:’** 

LINE 

20 MAY 
PIXEL 

1978 

ID 

HCMM 

VIS 

CCT DIGITAL NUfir-ERE 
DIR NIR DEL T 

OF POOR QUALlTf 
T I DIR/DEL 

POWER 

STATION 

COOLING PONDS MASON VALLEY. YERINGTON 

DISTRICT 

367 

682 

2 

POWR 

27. 0 

109 0 

35 0 

122 0 

54 0 

0. 893 

367 

683 

2 

POWR 

29. 0 

106 0 

36 0 

116 0 

56 0 

0 914 

367 

684 

2 

POWR 

39. 0 

119. 0 

31. 0 

136. 0 

43 0 

0. 875 

367 

688 

2 

POWR 

42. 0 

134 0 

26 0 

155. 0 

34 0 

0. 865 

368 

682 

2 

POWR 

26 0 

109. 0 

36 0 

121. 0 

55. 0 

0. 901 

368 

683 

2 

POWR 

29. 0 

no. 0 

41. 0 

116 0 

57. 0 

0. 948 

368 

684 

2 

POWR 

39. 0 

134 0 

33 0 

147. 0 

38. 0 

0. 912 

368 

685 

2 

POWR 

47 0 

146 0 

27 0 

165 0 

31 0 

0. 665 

369 

682 

2 

POWR 

33. 0 

116 0 

31 0 

133 0 

45 0 

0. 872 

369 

683 

2 

POWR 

33 0 

124 0 

32 0 

139. 0 

42 0 

0 892 

369 

684 

2 

POWR 

41. 0 

140. 0 

31. 0 

155. 0 

35 0 

0. 903 

369 

685 

2 

POWR 

47. 0 

147 0 

27. 0 

166 0 

31 0 

0. 866 

LEACH 

PONDS AND ORANODIORITE WASTE DUMPSr YERINGTON MINE 

386 

665 

1 

DUMP 

34. 0 

118 0 

40. 0 

124 0 

50 0 

0. 952 

386 

666 

1 

DUMP 

41. 0 

125. 0 

39. 0 

132 0 

44 0 

0. 947 

386 

667 

1 

DUMP 

44. 0 

125. 0 

40. 0 

131. 0 

44. 0 

0. 954 

386 

668 

1 

DUMP 

45. 0 

121. 0 

35. 0 

133. 0 

43. 0 

0. 910 

386 

669 

1 

DUMP 

42. 0 

123. 0 

32 0 

138. 0 

41. 0 

0. 891 

387 

665 

1 

DUMP 

28. 0 

107. 0 

42 0 

1 12. 0 

61. 0 

0. 955 

387 

666 

1 

DUMP 

39. 0 

109. 0 

43 0 

113 0 

57. 0 

0. 965 

1 387 

667 

1 

DUMP 

44. 0 

no. 0 

44. 0 

113. 0 

56 0 

0. 973 

1 387 

ib&E 

1 

DUMP 

44. 0 

1 16 0 

39. 0 

124 0 

46 0 

0 935 

387 

669 

1 

DUMP 

41. 0 

135. 0 

35. 0 

146. 0 

38 0 

0. 925 

388 

665 

1 

DUMP 

23 0 

101 0 

38 0 

111.0 

64. 0 

0 910 

388 

666 

1 

DUMP 

41 0 

97. 0 

40. 0 

105 0 

64. 0 

0. 924 

388 

667 

1 

DUMP 

56 0 

107 0 

40. 0 

114 0 

50 0 

0. 939 

3E8 

668 

1 

DUMP 

57. 0 

112. 0 

39. 0 

120. 0 

47. 0 

0. 933 

388 

669 

1 

DUMP 

49. 0 

127. 0 

37. 0 

136 0 

40 0 

0. 934 

389 

665 

1 

DUMP 

52. 0 

116. 0 

37. 0 

126. 0 

44. 0 

0. 921 

339 

666 

1 

DUMP 

57. 0 

114. 0 

39 0 

122 0 

45. 0 

0. 934 

389 

667 

1 

DUMP 

67. 0 

114. 0 

40. 0 

121. 0 

43 0 

0. 942 

389 

666 

1 

DUMP 

63 0 

118. 0 

38 0 

127. 0 

41. 0 

0. 929 

389 

669 

1 

DUMP 

52. 0 

126. 0 

37. 0 

137. 0 

39. 0 

0 934 

390 

665 

1 

DUMP 

50. 0 

120. 0 

38 0 

129. 0 

44. 0 

0. 930 

390 

666 

1 

DUMP 

54 0 

116. 0 

40 0 

123 0 

46 0 

0. 943 

390 

667 

1 

DUMP 

60 0 

118 0 

40 0 

124. 0 

43 C 

0. 952 

390 

666 

1 

DUMP 

60 0 

124. 0 

39. 0 

131 0 

40 0 

0. 947 

390 

669 

1 

DUMP 

47. 0 

135. 0 

36. 0 

144 0 

37. 0 

0 937 

391 

665 

1 

DUMP 

43. 0 

119. 0 

38 0 

126 0 

46 0 

0. 930 

391 

666 

1 

DUMP 

4). 0 

115. 0 

42. 0 

119. 0 

51 0 

0. 966 

391 

667 

1 

DUMP 

45. 0 

128. 0 

43. 0 

130. 0 

44 0 

0. 985 

391 

668 

1 

DUMP 

45. 0 

137. 0 

40. 0 

142. 0 

39. 0 

0. 965 

391 

669 

1 

DUMP 

45. 0 

135. 0 

36. 0 

144 0 

38. 0 

0. 937 

392 

665 

1 

DUMP 

41. 0 

127. 0 

38. 0 

135. 0 

43. 0 

0. 941 

392 

666 

1 

DUMP 

38. 0 

120. 0 

39 0 

127. 0 

47. 0 

0. 945 

. 392 

667 

1 

DUMP 

46 0 

130. 0 

40. 0 

135. 0 

41. 0 

0. 963 

i 392 

668 

1 

DUMP 

45. 0 

134. 0 

38. 0 

141. 0 

2®. 0 

0. 950 

392 

669 

1 

DUMP 

42. 0 

124. 0 

35. 0 

136 0 

42. 0 

0. 912 


HCMM REGISTERED TAPES RECEIVED TO NOV 1980 


ORDERED 

A-0034 31 May 1978 
A-0087 27 July 1978 
A- 0226 8 Dec 1978 


RECEIVED 


A-0024 20 May 1978 
(received Aug 1980) 
A-0082 22 July 1978 
(received Aug 1980) 
A-0226 8 Dec 1978 

(received Nov 10,1980) 


COMMENTS 

SITE WAS ON TAPE EDGE 
SITE NOT ON TAPE 
SITE UNDER CLOUD 


278.0 *K and correspond to areas surrounding the water bodies, the Yerington 
mine, etc. Class C are temperatures from 277.6 *K to 273*8 *K. AM the 
rest (cooler areas) are represented in white. Vegetation and high topo- 
graphic areas are included in Class C and in white. 


The following chart summarizes the quality of the data used in the 


study: 


August 

8 - 9 . 1977 


May 30, 

1978 


Oecembe r 

8 - 9 , 1978 


Field 

Measurements 


. Good . . • 

. • . . • . 


Black bodies 
Reversed 
Temp, problems 


Cal i brat ion 
Problems 


(Not Launched 
Yet) 


None taken 



Cal i brat ion 
Problems 


Poor due to 
weather conditions 
(freezing) 




Poor (no information 
contained in data) 


Analysis of the chart indicates that the field measurements corresponding 
to August 7-8, 1977 and the HCMM satellite data corresponding to May 30, 1978 


are the only data sets of use for the present study 






VI. CONCLUSIONS 


This study rslstss sn attsmpt to ostch ths HCMM and U2HCMR* 
darlvad tamparatura data ovar two taat altaa (ovar vary local slat) to 
similar data collactad in tha fiald at naarly tha same timas. Con- 
sidarabla loglatical problams vara ancountarad. Tha rasults indicate 
that HCMM invaatigations (using raaolution calls of 500 Oi or so) are 
bast conductad with araally-axtansiva aitas, rathar than point obaer- 
vationa. Tha DAY -VIS Imagary is of axcellant quality* and has con- 
siderable usefulness for GEOLOGY* especially for structural (lineament) 
studies. For these purposes one does not need the Day-Night registered 
imagery* except that as a single product* not to ba used for further 
calculations* the DELTA-T imagery is most useful* again for structural 
geology. Our attempts to register the ground-observed temps, (even for 
0.5 sq. mile targets) were unsuccessful* due to the excessive pixel- 
to-plxel noise on the HCMM data. 

Several computer models were explored* and related to changing of 
the values of thermal parameters* with observed data. Unleas quite 
complex models* with many parameters which can only be observed (perhaps 
not even measured!) under remote sensing conditions (e.g. roughness* 
wind sheer* etc.) the model outputs do not match the observed data. 
Empirical relationships may be most readily studied. 
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APPENDIX 1: DATA LISTS 

The following pages contain the data collected from the soil tempera- 
ture probes, sol-a-meter, net radiometer, PRT-A and PRT-5, air temperature 
sensors, psychrometert , anemometers, and cloud cover observations. 

The data is also labeled except f^r some of the PRT results. The top 
row for each time in the sandbox and spike data sheets show five voltage 
readings taken over a short period of time followed by the average of these 
voltages. The second row under column 1 is the temperature reading from 
the PRT-5 direct readout dial during the first measurement. The number 
under AVR6. in the second row is the temperature corresponding to the aver- 
age millivolt reading. 

PRT-k measurements were made at the MacArthur station, but not until 
2100, 8 August 1977 because of operational problems. All readings were 
recorded in degrees F. The second row under AVRG. is simply the average 
in *C. Several readings are listed as greater than ki.yc. This is a 
result of the surface temperature exceeding the range of the PRT-A. 
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APPENDIX 2: CALI lA^T ION OF THE SOL-«-METERS 

TabI* 2*1 fthows raadirigt takan from tKa Imo »ol *a'mttar» at Idantical 

timaa In naarly idantical condition*. Tha first tat of columns show tha 

rasults for matar I6S^07 Mhich Mas cailbratad by tha manufactura In 1966 

(Figura 2‘'1). Tha first column shows tha voltaga output of this sol*a* 

matar, usad with an ampllfiar that incraasad tha voltaga 660 timas. Column 

2 lists tha output In mV dividad by tha amplification. Tha naxt numbar 

across tha row is tha numbar of 9TU/Hr/SqFt takan diractly from llna 1 in 

2 

Figura 2-1. This valua is than convartad to cal/HW/cm using tha convar> 

tion factor shown » Tha last two columns for this 

unit ara basad on tha llna labalad 2 in Figura 2*1. 

Raadlngt for sol*a*matar 12102 ara Itstad in tha last thraa columns. 

Tha first of thasa columns shows the voltaga output without a.'i amplifier. 

2 

Using Figura 2*2, this numbar Is than convartad to cal/min/cm . This 
numbar is comparable to tha readings of tha same units for sol*a*matar 
I6S607. A quick comparison shows that line I In Figura 2*1 gives rasults 
more Ilka that of Figura 2*2 than does line 2 in Figura 2*1. Tha last 
column in Table 2*1 shows the difference between the two values. Beneath 
this numbar Is the purcant error in terms of tha 1966 values. 

Figure 2*3 is a graph of tha percent error at a function of tha sol- 

a*metar 165607 readings basad or Table 2*1 values. Tha graph serves as a 
way to correct the values of sol*a*metar 65607 In terms of tha more recently 
calibrated sol*a*matar 2102. 

Tabic 2*2 lists the original and corrected values from both sol*a~meters 

during the 26 hoir collection period. The last column gives the converted 

chart value, tSic corrected value obtained using Figure 2*3* Figures 8 and 9 
give these numbers plotted as a function of time. 
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APPENDIX 3: RECORDING THERNOHETERS 

Table 3*1 shows the readings for both recording thermofflcters A and 8 
under nearly Identical conditions. These readings are compared to those of 
the new 5810 thermomett unit. Asterisks indicate instances where the 
needle of recording thermometer A beceme stuck at higher temperatures. 

Averages of dl of these comparisons are also listed in this table. 

It is important to note that the probes were improperly set in the 
field. The whole probe and several inches of the wire should have been sub* 
merged Into the water because of the conductive metal wire connected to the 
probe. Instead I only >he probe was in the water. As a result, the recorded 
water temperatures a pear too variable due to the greater temperature change 
In the air. Due to this error, this attempt to calibrate the recorded data 
may not be very meaningful. 

Figure 3*1 is a plot of the temperature recorded by thermometers A and 8 
as a Function of the temperature given by the 5810 unit. Lines were visually 
fitted to the points. Note that the two circled points represent values found 
in the field. Equations of these lines were then celculnted for both thermome* 
ters. 

Figure 3*2 is sieply a graph of temperature readings from thermometers A 
versus B. The points represent values recorded for both thermometers in the 
same water bath. 

Table 3*2, 3*3, end 3*^ ere the temperature recording charts made at 
Anaconda's leach pond. 

Table 3*5 lists the temperature values for both the thermometers every 
30 minutes. The correction factors are from Figure 3*1. The corrected 
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temperatures ere also l.sted for each original vAlue of A and 6. Figures 1i 
and 12 show these values plotted as a function of time. The bottom of T^ible 
3*5 lists the temperature range of hoth stations. The first calculation is 
based upon raw data from the recording chart, the second set is based on the 
correction terms used in this report. The final listing uses the calibration 
factor determined in the field (circled points in Figure 3*1). 


TABLE 3*1: Calibration of Recording THenmetert 
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APPENDIX k: SOIL SAMPLES 

Table k~\ lists the wet and dry welg. ts of each soil samp'e, and the 
percent of moisture by weight. Each sample is labeled by station (Anaconda 
or MacArthur) and by depth below the surfac*. Asterisks indicate those 
samples which were determined to be too small tc be employed. Plate 5 shows 
samples after the drying process. 

Samples were collected at two different times. One set was obtained 
August 1977 while the hole was being dug out. Th( soil in the auger ' is 
deposited into a plastic air tight sample bottle. Additional samples were 
collected August 9, 1977 after thermal measurements were completed. These 
samples were obtained while dismantling the st». ^ns. 

Each sample was weighed before and after being dried in an oven. Samples 
were removed from the sample bottles and placed on previously weighed watch- 
glasses. Unfortunately, the very dry soils were gaining weight very quickly 
dui ^ig the time taken to weigh the samples. Thousandths of a grar change were 
rc.dily discernible for both dry samples gaining weight, and moister soils 
losing weight. In each case, those which were losing weight had a higher 
water content than those which were gaining weight. No weight change was 
noticed for the smaller samples collected August 9, 1977 from the MacArthur 
sewer pipe. The wetter samples were re-weighed after airdrying for A5 minutes. 
Changes of over 0.1 gram were observed. 

Each sample was very carefully place*' 'nto a drying oven, made by the 
National Appliance Company. The oven was set at approximately A.S which 
corresponds to a temperature range of 112- 120*C. All samples were allowed 
to dry over 2A hours. To ascertain the dryness of the soil, one can weigh 


no 


the samples periodically during the drying period until no discernible change 
in weight is observed. Warmer temperatures and more time can aiwuys result 
in a drier soil, if only by breaking down some of the clay minerals in the 
sol I . 


TABLE A-1: Moisture Content of Neveda •'ThenMl Test** Solis Collected August 5- 9, 1977 
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Plat* Sc. soa SAMPLES AFTER OVER DRYING. MOTE THE COARSE TEXTURE 
OF THE SAMPLE IM LEFT CORNER, FROM THE DUMP SITE SURFACE 
BY THE "SBWER”PIPE INSTALLATION. NOTE ALSO THE FINE YELLOW SAMPLE 
JUST BENEATH THE COARSE SAMPLE. THIS SOIL HAS FROM THE DUMP SITE 
SEWER PIPE LOCALITY 4.5 INCHES BELOW THE SURFACE, AND WETTER THAN 
THE OTHERS. 
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APPENDIX 5: CALIBRATION OF THE PRT-1* AND PRT-5 UNITS 

To calibrat* the PRT instruments a black body filled with well mixed 
water was used as the sensing target. The temperature of the water measured 
by a reliable thermometer probe was compared to the PRT reading of the water 
temperature. 

Table S'l lists temperatures taken v/ith this equipment using a mercury 
thermometer to measure the water temperature. Voltmeter readings from the 
PRT are recorded in the column following the PRT'-S meter readings. The 
letter after the voltage indicates the srcle of the reading; high, medium, or 
low. The next column lists the temperature equivalent of the voltage. The 
difference between the temperature from the PRT-5 meter and the equivalent 
voltage temperature is listed in the next column. The last two columns show 
the difference between the PRT-5 meter readings and the voltmeter temperature 
readings. 

Figure 5*1 is a graph showing the difference petween the voltage tempera 
tore of the water and the measured water temperature as a function of the 
voltmeter temperatures. The two curves show the trends for the medium and 
high range. 

Figure 5*2 is a plot of the PRT-5 voltage temperature versus the water- 
bath temperature foV both the medium and high range readings. 

An attempt was made to determine the effects of the ambient air tempera- 
ture by comparing readings with hot ali blasting the PRT-5 sensor head with 
cooler air surroundings. No effect was detected. While running this experi- 
ment, a new set of data was recorded. Table 5*2 summarizes this new data 
set in the same format as Table 5*1. Figures 5*3 and 5*1* present this new 
data. 
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The PRT-A was calibrated the same way as the PRT-5. Table 5"3 summarizes 
the resulting data. However, no tables were available to convert voltages 
into temperature. Itote also that readings within 2*F were as accurate as could 
be made. 

Figure 5*5 Is a plot of the PRT-A voltmeter readings as a function of the 
temperature readings of the waterbath. Figure 5*6 attempts to relate PRT'A 
meter readings to the voltage readings. The relationship is described by the 
equation of the plotted line. 

The accuracy of the PRT-I* may be shown by Figure 5*7 which compares the 
PRT*4 readings to the measured water temperature. Table 5*^ lists the call* 
bration data obtained the evening before the field study was begun. The 
format of this chart is the same es for the previous tables. 
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TABLE 5*2: Calibration of PRT-5 Meter Readings 


Meter 

•C 


Volt- 

Meter 

V 

4.987h 

♦ . 848 h 
♦.828h 
♦.793H 
♦.705h 
+.670h 
♦.574h 
♦.54lh 
4,i«44h 
♦.360h 
♦.321h 
4.26lh 
4.238h 
+ .213h 

♦ .173b 

♦ .152h 

132h 

♦ .lOlh 

♦ .094h 
+ .924m 

♦ .056h 
♦ . 886m 
+ .032h 

898 m 
+ .000h 
♦ . 809m 
■f . 804m 

♦ .777m 

♦ ,748m 
♦. 728 m 
♦,644m 

♦ ,651m 

♦, 629 m 
♦, 606 m 
♦.577m 
♦. 558m 
♦. 538m 
♦,503m 
♦,482m 
♦,469m 

♦ ,445m 
♦, 430 m 
+ ,402m 
♦,370m 

♦, 347 m 

♦. 333 m 

♦,314m 


Volt- 

Heter Voltmeter 

•C 


74.6 

69.7 

66,9 

67.7 

64.8 

63.6 

60.4 

59.2 

56.1 

53.1 

51.8 

49.6 

48.8 

47.9 

46.5 

45.8 

m 5.0 

43.9 

43.6 

42.4 

42.2 

41.1 

41.7 

40. 1 

40.0 

38.4 

38.3 

37.4 

36.4 

35.8 

34.7 

33.2 

32.5 

M.7 

30.7 

30.1 

27.4 

28.3 

27.6 

27.1 

26.3 

25.8 

24.8 

23.7 

22.8 

22.4 

21.8 


0.2 

♦ 0.8 

♦0.9 

♦ 0.8 

♦ 0.8 

♦0.9 

♦ 0.3 
♦ 1.0 
♦0.9 

♦ 0.7 

♦0.7 

♦ 0.8 

♦ 0.7 

♦ 0.6 

♦ 0.5 

♦ 0.7 

♦ 0.5 

♦ 0.6 

♦ 0.6 

♦ 0.1 

♦0.6 

♦ 0.4 

♦0.3 

♦ 1.4 

♦0.5 

♦0.6 

♦ 0.5 

♦0.6 

♦ 1.0 
♦0.7 
♦ 0.8 
♦0.9 
♦0.9 
♦ 0.8 
♦ 0.8 
♦0.9 
♦0.9 

♦ 0.7 

♦ 0.8 

♦0.9 

♦0.8 

♦ 0.7 

♦ 0.7 

♦ 0.8 

♦ 0.8 

♦ 0.7 

♦ 0.7 


Thermometer 

(5810-1)*C 

74.2 

69.6 

68.6 

67.0 

66.1 

65.1 

62.1 

60.4 

57.1 

54.0 

52.4 

50.5 

49.5 

47.7 

46.1 

46.0 

45.1 

43.8 

43.6 

43.6 

42.0 

42.0 

41.0 

41.0 

39.6 

39.8 

39.0 

38.0 

37.1 

36.2 

35.3 

33.9 

32.8 

31.9 

30.0 

30.1 

28.9 

27.9 

27.2 

26.7 

25.7 
25.1 
24.0 

22.7 

21.8 

21.3 

20.5 


Thermometer 

*C 

♦ 0.6 
♦ 0.5 
♦ 1.2 

♦ 1.5 

- 0.5 
- 0.6 
-0.9 
- 0.2 
- 0.1 
- 0.2 
♦ 0. 1 
- 0.1 
0.0 
♦ 0.8 
♦0.9 

♦ 0.5 

♦ 0.4 

♦0.7 

♦ 0.6 

- 1.1 

♦ 0.8 

- 0.5 

♦ 1.0 
♦0.5 
♦0.9 
- 0.8 
- 0.2 
♦ 0.0 
♦0.3 
♦0.3 
♦ 0.2 
♦ 0.2 
♦ C .6 
♦ 0.6 
♦0.7 
♦0.9 

♦ 0.4 
♦ 1.1 
♦ 1.2 
♦1.3 

♦ 1.4 

♦ 1.4 

♦ 1.5 

♦ 1.8 
♦ 1 .8 
♦ 1.8 
♦ 2.0 
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TABLE 5“3: Calibration of PRT-A Meter Readings 


Date 

— — 

Meter 

of 

Vol t- 
Meier 
V 

Thermometer 

(5810-1)*C 

— — 
Meter 

•c 

Thermometer 

•c 

8/27/77 

> no 

•f.016 

56.3 

> 43.3 

M ^ W 

in lab 

> no 

♦ .015 

54.8 

> 43.3 

- - - - 

(A-16) 

> no 

+ .018 

53.1 

> 43.3 

m ^ m m 


> no 


52.6 

> 4.33 

m m ^ ^ 


> no 

♦ .028 

51.8 

> 43.3 

m m ^ m 


> no 

♦ .041 

48.9 

> 43.3 

— 


> no 

♦ .043 

47.8 

> 43.3 

m m ^ m 


> no 

♦ .051 

47.0 

> 43.3 

mim mt m 


> no 

♦ .054 

46.1 

> 43.3 

m m m 


> no 

♦ .064 

45.0 

> 43.3 




> no 

♦ .071 

43.4 

> 43.3 

m m m 


> no 

♦ .080 

41.8 

> 43.3 

mmmm 


107-109 

♦ .093 

39.2 

42.2 

♦3.0 


lOA-106 

♦ .096 

38.2 

40.6 

♦2.4 


103-105 

♦ .100 

37.4 

40.0 

♦2.6 


102-10A 

♦ .104 

36.4 

39.4 

♦3.0 


100-102 

♦ . 109 

35.4 

38.3 

♦2.9 


97-99 

♦ .114 

34.0 

36.7 

♦2.7 


96-90 

♦ .119 

32.9 

36.1 

♦3.2 


9A-96 

♦ .122 

31.9 

35.0 

+3.1 


92-9^ 

♦ .128 

30.8 

33.9 

♦3.1 


91-93 

♦ .132 

30.1 

33.3 

♦3.2 


90-92 

♦ .135 

28.9 

32.8 

♦3.9 


88-90 

♦ .140 

27.9 

31.7 

♦3.8 


86-88 

♦ .142 

27.1 

30.6 

♦3.5 


86-88 

♦ .144 

26.7 

30.6 

♦3.9 


8i»-86 

♦ .150 

25.7 

29.4 

♦3.7 


8A-86 

♦ .151 

25.1 

29.4 

♦4.3 


82-8A 

♦.154 

24.0 

28. 3 

♦4.3 


80-82 

♦ .157 

23.4 

27.2 

♦3.8 


79-81 

♦ .161 

22.6 

26.7 

♦4.1 


78-80 

♦ .163 

21.8 

26.1 

♦4.3 


77-79 

♦ .166 

21.4 

25.6 

♦4.2 


75-77 

♦.168 

20.5 

24.4 

♦3.9 


7A-76 

♦ .170 

19.8 

23.9 

♦4.1 


7^-76 

♦ .175 

19.1 

23.9 

♦4.8 


73-75 

♦.173 

16.8 

23.3 

♦4.5 


72-7A 

♦ .173 

18.4 

22.8 

♦4.4 


70-72 

♦.179 

17.2 

21.7 

♦4.5 


70-72 

♦ .180 

16.6 

21.7 

♦5.1 


68-70 

♦ .183 

16.1 

20.6 

♦4.5 


68-70 

♦.184 

15.6 

20.6 

♦5.0 


67-69 

♦ .187 

15.1 

20.0 

♦4.9 


66-68 

♦ .187 

14.5 

19.4 

♦4.7 


65-67 

♦ .190 

14.0 

18.9 

♦4.9 


6A-66 

♦ .192 

13.6 

18.3 

♦4.7 


62-64 

♦ .196 

12.7 

17.2 

♦4.5 


62-64 

♦ .196 

12.2 

17.2 

♦5.0 


89-90 


mi rror 

32 










CAulBR/^TION PRT-4 VOLTMETEP^ 
Tei7\PBPATURB RgAO»NaS 
8/27/77 
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TABLE 5~*»: Calibration of PRT-4 and PRT-5 
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APPENDIX 6 

PHOTOGRAPHS OF FIELD WORK AUGUST 1977 



ORIGINAL PAGE 

BLACK AND WHITE PHOTOGRAPH 



Plat* 5a. LOOKING S OF LEACH POND WITH HECORDING THERMOMETERS, TOWARDS A 
HEMATICALLY-STAINED TAILINGS PILE 



Plat* 5b. WALKER RIVER NEAR YERINGTON, NEVADA 


i 


ORIGINAL PAGE 

BLACK AND WHITE PHOTOQMfM 
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REPORT ON THERHAL-INERTIA RELATED COMPUTER PROGRAMS 
FS6 6 FSFL2 (CSIROTEMP) 

WATSON, SURTEMP 

by 

Randy Kirk 

Stanford Ramota Sansing Laboratory 
Stanford Uni varsity 
Stanford, California 91«30S 

It If difficult to compara tha four programs on the basis of eomplaxity, 
for they differ widely In purpose: FS6 attempts to fit observed temperature 

data by adjusting two or three parameters (with others held fixed), while 
WATSON merely calculates diurnal temperature curves corresponding to a series 
of Input parameters, any given one of which may be assigned multiple values for 
a family of curvas. The machinery required to allow selection of any parameter 
as the 'variable' with multiple values makes WATSON appear much more complex 
than It otherwise would. SURTEHP and FSFL2 marely accept a full set of Input 
parameters and calculate a temperature curve ; they lack as well the plotting 
routines shared by the other two programs. FS6 reads data from and writes to 
text files, while the other programs deal interactively with the. keyboard 
user. 

FS6, FSFL2, and WATSON employ modifications of the same algorithm, fourier 
decomposition of the solar heat flux, linearization of other heat flux terms, 
and fourier representation of the surface temperature. They differ In that 
FS6 attempts to fit obeserved temperature data whereas WATSON merely calculates 
one or more diurnal curves, and as well in the complexity of their parameter 
sets. WATSON Includes the effects of emissivity, cloud cover, and geothermal 
flux, while FS6 lumps all nonsolar flux terms into a single variable. 

FSFL2 is an earlier, nonfitttng version of FSS, which includes the effects of 
sensible heat transfer and longwave radiation from both earth and sky, but not 
cloud cover or slope. SURTEHP has a parameter set similar to WATSON, but 
utilizes an entirely different algorithm. The Laplace transform method of 
Jaeger is used to express the total thermal flux required to maintain a 
periodically recurring temperature curve (discretized into a twenty-point series) 
on the basis of heat conduction In the earth. Models of soiar flux, longwave 
emission, and so on give a second expression for the fiuxes as a function of the 
temperatures at the twenty times. A least-squares fitting routine adjusts the 
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t«mp«raturt until the two expersslons are consistent. A set of observed 
Temperatures may ba read In and printed out for comparison, but are not 
used In the calculation. It is worthwhile to emphaslie that, thouph FS6 
end SUATEHO Util lie the same I east “squares routines, they use them for entirely 
different purposes. In SURTEHP, they find a set of temperatures which satisfy 
the heat flux requirements for periodicity, given certain other parameters such 
as thermal Inertia. In FSb, they operate on the temperature curve as a given 
vector-valued function of the thermal inertia and other parameters, which they 
adjust to modify the solution. 

A brief description of the computations Involved In each program Is given 
below, followed by a discussion of the various mathematical approaches to the 
heat equation problem. I have attempted to convert the symbology of the Jaeger 
and Watson papers and the three programs to the following consistent set: 


Thermal Inertia 

Albedo 

Emissivlty 

Air Temperature 

Sky temperature (mean) 

Sky temperature (day) 

Sky temperature (night) 

Latitude 

Solar Declfnwtior 
Solar zenith angle 
Local surface dip 
Strike 

Atmospheric transmissivity 

Air pressure 

Water content 

Dust oiitent 

Von Kerman's constant 

Wind Friction velocity 

Surface roughness length 

Height of meterologic observations 


P 

A 


c 



X 

6 

Z 

d 


s 

H 

P 

W 

CD 

k 



133 


StabM j ty-depend«nt profile function 
Heen turfect'elr tompereturc difference 
ftetio of tky to earth emlttlvity 
Free eky fraction 
Solar constant 
Observed ground tamparatura 
Calculated ground tamparaturas 
Diurnal angular frequency 
Time (from noon) 

Total heat flux 
Geothermal flux 
Solar flux 

Fourier amplitude t phase of solar flux 
Sine and cosine amplitudes of solar flux 
Thickness of layer over the half-space 
Thermal inertia of covering layer 
DIffusIvity of covering layer 


if 


c 

C 

T 

V 


b> 

t 

F 

Q 

I 



P 


I 


K 
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Calculation proeaadt axactly as in fS6, axcept that tha llntar heat'transfar 
tarm includes not only upward longwave radiation, but downward radiation 
and sensibla heat loss as well, using the assumptions that air temperature 
is a linear function of ground temperature, and sky temperature proportional 
to air tamparatura: 


6 - W(l-B) ♦ CO (kV^ - (V-AT)^) 


0- 


In ( 0.8905 

ln(0.8905 i4«^/kU^) ) 


so that 

Once the values of temperature V are known, rather than adjusting them by 
altering P, G, and V as in FS6, the program calculates and prints each 
of the terms in the energy budget for each point In time: 

\lr ■ ^alr 8(V(t)-V), l(t), l(t) ♦ co(ctJ,^ - V^), 

H(t) - - V), and F(t). 


WATSON 


(latest version Implemented by T.E. Townsend) 


Input parameters: P, A, e. A, 8, *» Q, C 

Any one parameter may have 2-5 values, while all others have one. 


The solar Irradlance is modeled by the equations: 
M* 1 - 1/5* (cos(A)cos(8)cos(u>t) - sin (X)sln(6))^ 


cos(Z, ■ cos(X)cos(6)cos(d)cos(b>t) ♦ cos(X)sln(6)sln(d)sln(<Jt) 


- sin(X)cos(5)sin(d)cos(s)cos(:i)t) - sin(X)sin(6)cos(d)cos(s) 


I ■ Mcos(z) 


and avaluatad at 99 aqually-»paced points In timo to ylald by racursiva 
calculation a SO~hartnonlc pourlar cos I na>and-s I na sarlas: 
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$9 

l(t) • y cos(rvt)t) ♦ B_ sln(nu>t) 
o ••jk n n 

n*0 


l’ 

n«0 




inu)t 


Tha Fourier coafficlants for tha tamparatura ara calculatad by tha complax*’ 
admtttanca mathod of Byrne and Davls» parmittinq Inclusion In tha modal of 
a layar ovar tha half-spaea with diffarant tharmal propart las. Tha calcula- 
tion Is most simply axprassad using tha complax-axponentlal raprasantatlon of 
tha Fourlar sarlas: 


V(t) 



wher. 


— V - 

♦ F a 


♦ (I -A) S^C 
Iit/I» 



♦ '»n> 

7 

n 



^Inut 


in tha casa of no layer, 
P| « P, I • 0, and 


Tha maan tamparatura Is calculatad by the expression: 
V - (l-A)A C ^ ♦ f ♦ 


sky 


(1'A)S CA^ ♦ Q 
o o 

3 — 


♦ T 


Than for l»9 equally-spaced times the temperature is evaluated: 

<[9 


V(t) 


kcofj 


sky 


® n *"* (^"”1 (^ b ) t ) 
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♦ <»„«• 



-A l« 

n 



t i n (nai* ) ) 


•nd tha r«»ultlng t«mptr«tur« curve plottr,d; the procest It repeated for 
•II veluet of the varying parameter. 
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SURTEMP 

(oldest vrtton of WATSON programs) 

Input p'irameters: P, A, C, Tp, c, X, 6, d, s, (t, T(t)) pairs 

The values of T will be printed out with the result V, but ere never 
Involved In calculation. 

The model calculates all values at an evenly spaced set of twenty 
times, starting with the solar flux: 

M ■ I - 1/5 (cos(X) cos(6) cos («t) + sln(X) sln(6))^ 

cos(X- 6 sln(s)) cos(6) co$(a)i ♦ d cos(s)) sln(X-6sln($))sln(6) 
* " or zero If this Is negative 

T . ^I>, cos(Z)> 0 
T^. cos(Z)< 0 

l(t) - ^ C(I-A)M cos(Z) + CoT^Jy 


The Jaeger method gives a prediction of the total heat flux Into the ground 
required to sustain a periodic temperature curve: 

p 20 

f, y V, where Y, ere constants 

' ^ * Ji1 • ' ‘ 

The program defines an error In the flux estimation: 


AF 


I 


♦ oeV 


I 


end this Is fit to zero as a function of the Vj's by a least squares method. 
The data re yped and pseudoplotted on the terminal. No options to fit exist, 
but the parameters may be adjusted and the calculations repeated. 


FOURIER SERIES SOLUTIONS TO THE HEAT EQUATION: 
THERMAL INERTIA. THERMAL ADMITTANCE. AND THE C-FACT OR 
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Temperature -curve analysis Is concerned with the solution of the heat- 
conduction equation 



d „ 3T 

3* at 


at 


0 ) 


In the earth, with B.C. at the upper surface 

I* - ♦ 6 ♦ H LE - 0 (2) 

tot 

where I. ^ total radiative flux, G > heat conduction Into the earth, I.e. 
tot 

the heat flux In (1), H ■ sensible heat loss to the air, and LE ■ evaporative 
heat loss. 

Solutions for the simplified case H * LE * 0, I. . ■ A cos(u)t) * B sln(u)t) 
2 2 4 “IB 

■(A -f B )\os(wt - tan ^ ) ■ (A-IB) exp (l(wt)} Is possible; we get 


2-0 


. A - IB 


(D P 


exp { I (o)t " ij-) } 


where P ■ C Is the thermal Inertia. Dividing the temperature Into 
the flux we can obtain the admittance (by analogy with admittance In A.C. 
electronics, where current corresponds to flux, and voltage to temperature). 


Y - 




(3) 


Note that Y Is a function of both w and P, and that It Is proportional 

to e^ since T lags I by The In-phase and quadrature components 

of 1 are equal In magnitude. 

in the full problem, I Is Just one of the terms in the boundary condition. 
*to'i. "driving" sinusoidal flux, but upward thermal 

radiation and downward sky ratlatlon. Sensible heat transfer and evaporation 
also operate in parallel with radiation. Unfortunately, though all of these 
depend on temperature, none of them express a simple (but complex. In the 
sense of Involving I) proportionality: 


PRECEDING PAGE BLANI^ NOT FILMED 
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'tot ♦ H ♦ LE - YT ♦ l^(T) ♦ (-I^(T)) ♦ H(T) + LE(T) 

The solution Is to linearize the dependencies of these quantities, I.e. to 


assume for example that l^(T) - 


T and so on. The constant term Is 


To 


removed by virtue of the linearity of the differential equatlon--lt merely 
adds a constant amount to the mean temperature. Thus, we define an effective 
ad'fil ttar.ce : 


r 


H + LE 




TY 


eff 


Since these functions depend on the Instantaneous temperature they are real and 
Independent on u). They are: 


lit 

3T 


ili 

3T 


3H 

Tf 




lie -? 

^ sky° sky 3 T 


U^kpK 


ln(z/z^-AV>j^) 


3^1 r . 
3T ^ 


(5a) 

(5b) 

(5c) 


3E 

Tt 


3H 

5T 



3S 



(5d) 


(k ■ vonKarman's constant, • friction velocity ■ roughness length, 

2 ■ height of evaluation of ■ stability parameters, S ■ humidity). 

The *G factor* Is the sum of these effective admittances, or some of 
them. Various authors use different symbols for it and include a different 
selection of terms. In any case, one gets 


Vf " 6-^Y-G+vCpeJ-(G +\/| P) ♦ i (V| P) 


( 6 ) 
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^ ^ A ■ IB ^Iwt 

(G 4 /f P) ♦ l(7y P) 


( 7 ») 




•*• (-a/| P - (G +yf p) B) sin (wt)} 


(7b) 



cos(ut ' tan 



-tan 


-1 


f V 

I -T^TJfrrfT-J ) 


in the three types of expression of the Fourier term. The solution to 
the heat equation for a realistic diurnal temperature curve may now be 
assembled by decomposing the solar flux Into a series of harmonic terms 
and utilizing equation (7) to find the temperature variation due to each 
component. 


141 


Other Approaches to the Problem 

Slmptlf teetton of the boundary value problem for the heat equation 
to the point where an analytic solution In the form of Fourier series 
Is feasible Is not the only method In the literature. Three other 
approaches have been used. 

Jaeger utilized Laplace transform methods to express the total heat flux 
out of a semi 'Infinite body (the ground) required for a given temperature 
profile to repeat Itself periodically. He divided the period •* the "day" - 
Into a number of discrete Intervals* and found the flux In the 1^^ period 
^1 for a temperature of unity In the first period and zero in the others. 

Then a complete temperature curve could be built up by summing over a set 
of these temperature "pulses" of different strengths at different times: 

( 8 ) 

This is an expression for the flux In the 1*^ Interval In terms of the 
temperatures in the other intervals Tj. The computation rests on the assump- 
tion that the fluxes for different temperature "pulses" may be summed, that 
is, that the problem Is I inear , the same assumption made in the fourier 
solution when harmonic terms of different frequencies were summed. In 
fact, the boundary conditions are nonlinear, but can be approximated by 
linear functions over the small temperature range within each time Interval. 
The full boundary conditions yield a second expression for the fluxes 
in terms of the temperatures: 

F| - F(T,) - I, + l| (T,) + (-1^ (T,)) + H(T,) ♦ Lr.(T,) (9) 

Where, as above, the terms on the right are Insolation, upward and downward 
longwave radiation, sensible heat transfer, and evaporative 'teat transfer. 

The thermal curve is computed by equating (8) and (9) for all Fj, and seeking 
the T| for which the equality holds by some approximate methods. 

A very different method which has been used by Rosema and Fratt, among 
others. Is to attempt to solve the partial differential equation (I) with 
BC (2) by purely numerical means, using some particular differential equation 
solving algorithm. In general, the differential equation is converted to 
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• difference equation by evaiuating the temperature on a discrete iattice 
of times and depths. An initiai temperature profiie is specif ied» and 
an Implicit or explicit finite difference relationship is used to generate 
successive profiles at each of the time intervals. The values at zero 
depth constitute the results of interest. 

A full discussion of the numerical solution of the heat equation 

Mould be out of place here« but some of the advantages and disadvantages of 

the method may be outlined. Mott significantly, any type of boundary 

conditions, linear, nonlinear, arbitrary, or a mixture of modeled and 

observed fluxes may be used with equal facility, similarly the soil proper* 

ties need not be assumed homogeneous; an arbitrary profiie may be modeled. 

As a final advantage, data are available in depth as well as at the 

surface. On the other hand, the choice of an algorithm, location of 

points in the vertical profile, and the choice of a time increment, 

represent a highly complex tradeoff between numerical stability, roundoff 

errors, and computation time. At best, numerical methods require far 

more computation than analytic ones. Furthermore, they yield only thermal 

curves for particular sets of physical parameters, whereas the fourier method, 

for example, yields a closed expression for the maximum variation in 

temperature as a function of thermal Inertia and so on. The analytic 

3T 

methods deal with an Implicit lower boundary condition '^'^0 as z , 

whereas an assumed temperature or geothermal flux must be specified at a 
finite depth in the numerical methods. An initiai condition must also 
be explicitly specified. In practice, an estimate of the initial profile 
is made and the heat equation solved over a period of one day to obtain a 
better, "relaxed" estimate with which the computation of interest are made. 

A final approximate method has been discussed by Deardorff. The 
heat equation (1) is split Into two first order equations: 

pAzC - - (f + (10a) 

6 - -K (10b) 

where F and are respectively the non-conductive and conductive fluxes 

out of a small volume of material. F vanishes except at the surface, and 
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the uifference between the downwerd-d tree ted fluxes G at the 

lOl 

bottom and top of the volume. The thickness of material eonstdared Is 
I. Now, F It modeled analytically, as we have te<rn above. If 6^^^ at 
the turface could be modeled alto, the tacond order partial differential 
aquatlon(l) would be converted to a firtt order ordinary differential 
aquation (lOa) and tolutlon would be relatively trivial. Oeardorff ditcusses 
a number of tuch modelt, finding that of Bhumralkar and Blackadar to be 
the mott effective. With Vr representing the ground turface temperature 
this method yields: 


^ - e, f/(pCd) - cj (11) 

Here ?2 Is the toil temperature at tome depth d where It remains fixed on 
the time scale of a few days (a similar calculation may be done to find 
te effect of the annual temperature variation). C, and C. are constants; 

4 ' ^ 

Bhumralkar sets ■ 2ir and C 2 ■ 2ir while Blackadar uses the vaiues 
3.72 and 7 . A respectively. Typical RMS errors for this method with respect 
to a full f ini te-difference tolutlon of (1) were found to be close to 
k% over a wide variety of soil and atmospheric conditions. The simplification 
involved in this method may thus be warranted by the great reduction in 
computation time with little loss of accuracy. 
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